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ABSTRACT 

Context. The large number and all-sky distribution of quasars from different surveys, along with their presence in large, deep astro- 
metric catalogs, enables the building of an optical materialization of the ICRS following its defining principles. Namely, kinematically 
non-rotating with respect to the ensemble of distant extragalactic objects; aligned with the mean equator and dynamical equinox of 
J2000; and realized by a list of adopted coordinates of extragalatic sources. 

Aims. The LQRF was built with the care of avoiding incorrect matches of its constituents quasars, homogenizing the astrometry from 
the different catalogs and lists in which the constituent quasars are identified, and attaining global alignment with the ICRF, as well 
as fully consistent source positions 

Methods. Starting from the updated and presumably complete LQAC list of QSOs, the initial optical positions of those quasars are 
found in the USNO B1.0 and GSC2.3 catalogs, and from the SDSS DR5. The initial positions are next placed onto UCAC2-based 
reference frames, following by an alignment with the ICRF, to which were added the most precise sources from the VLBA and VLA 
calibrator lists - when reliable optical counterparts exist. Finally, the LQRF axes are inspected through spherical harmonics, to define 
right ascension, declination and magnitude terms. 

Results. The LQRF contains 100,165 quasars, well represented across the whole sky, with 10 arcmin being the average distance 
between adjacent elements. The global alignment with the ICRF is 1.5 mas, and the individual position accuracies are represented by 
a Poisson distribution peaking at 139 mas in right ascension and 130 mas in declination. 

Conclusions. The LQRF contains J2000 referred equatorial coordinates, and is complemented by redshift and photometry informa- 
tion from the LQAC. It is designed to be an astrometric frame, but it is also the basis for the GAIA mission initial quasars' list, and 
can be used as a test bench for quasars' space distribution and luminosity function studies. 
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1. Introduction 

The establishment of the ICRS as the celestial reference system 
(Arias et al., 1995, Feissel & Mignard, 1998) answers not only 
the desiderata in terms of the choice of the physical model and 
coherence, but also of practicality when confronted with the ear- 
lier dynamical model. This is especially true since the system is 
materialized by the direction of quasi-inertial grid points, which 
are given by distant, extragalactic sources, instead of nearby, 
moving stars. However, for a long time after the establishment 
of the ICRS, the faintness of the extragalactic sources prevented 
a precise determination of optical positions for a large num- 
ber of them. The VLBI radio technique answered this need, by 
providing precise positions. But again, such precise positions 
could only be derived for a limited number of strong and sta- 
ble sources, and the ones with richer observational history have 
come to define the ICRF (Ma et al. 1998). Yet even this small 
number of sources enabled the VLBI to make important contri- 
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butions in defining both geophysics and radio astronomy stan- 
dards. 

In contrast, concerning extragalactic sources on the optical 
domain, poor astrometry and small number of known sources 
made with which the standard representation of the ICRF were 
assigned to the stellar catalog Hipparcos (Perryman et al., 1997). 
Currently several extensions of the HCRF (Hipparcos Catalog 
Reference Frame), densify this representation, reaching fainter 
magnitudes. 

The situation concerning the quantity of optically recognized 
bright quasars as well as concerning their sky distribution has 
witnessed significant changes in recent years (Veron-Cetty & 
Veron, 2006, Souchay et al., 2008). At the same time, the sci- 
ence's demand for accurate astrometry of these objects has also 
increased. Examples include micro and macro lensing, binaries, 
and space density counts, as well as the requirements of space 
astronomy missions (Andrei et al., 2008) and the very improve- 
ment of the stellar catalogs, that requires for a dense mesh of 
fiducial points (Fienga & Andrei, 2002, 2004). 
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Therefore, the conditions are put together to combine the 
large number of quasars optically recognized, with the large, pre- 
cise optical catalogs in order to produce an optical materializa- 
tion of the ICRS in terms of its first principles, namely kinemat- 
ically non-rotating with respect to the ensemble of distant extra- 
galactic objects, aligned with the mean equator and dynamical 
equinox of J2000, and realized by a list of adopted coordinates 
of extragalatic sources. This article presents one such realization, 
hereafter termed Large Quasar Reference Frame (LQRF). 

In the following section the input data is presented and re- 
viewed, whereas in its final subsection a schematic of the data 
stream is presented. The data stream is detailed in Sects. 3, 4 
and 5, where the analytical expressions used to derive the cata- 
log's intermediary steps are shown. In Sect. 3 a local astrometric 
solution is employed to place all quasar input catalogs onto the 
UCAC reference frame. In Sect. 4 the positions derived in this 
way are aligned with the ICRF frame. In Sect. 5 local inhomo- 
geneities are tackled by spherical harmonics fitting and local av- 
erage redressing. The LQRF catalog is presented in Sect. 6, and 
a summary and follow-up perspectives are presented in the last 
section. 



2. Input data and data flow 

The LQRF is built from gathering, comparing, and adjusting the 
positions of quasars and their stellar neighborhood, extracted 
from several catalogs and lists. A brief summary of their con- 
tent is presented in the following subsections. The role of each 
one is also mentioned, for the sake of introducing the general 
flow of the data treatment. The final subsection outlines the data 
flow, leaving the details for the following main sections. 



2.1. LQAC 

The Large Quasar Astrometric catalog (Souchay et al., 2008), 
hereafter LQAC, compiles the detections of 113,653 quasars 
from the original entries of the eight largest lists, with the sur- 
plus extracted from the latest edition of the Veron & Veron cat- 
alog (Veron-Cetty & Veron, 2006) after an analysis on the reli- 
ability of the input contributors. The LQAC contains extensive 
information on the physical properties of each object. In particu- 
lar, it contains, whenever available, the redshift, the magnitudes 
of the objects at 9 different optical bands (u, b, v, g, r, i, z, J, 
and K), and the radio fluxes at 5 different frequencies (lAGhz, 
23Ghz, 5.0Ghz, SAGhz, and 24Gte). The position of each ob- 
ject is given by the most precise equatorial coordinate directly 
available from the input lists. Therefore, the precision is quite 
variable, ranging from the sub-milli-arcsec of VLBI determina- 
tions to the level of a tenth of arcsec of the optical determina- 
tions, up to reach nearly one arcsec of single dish determina- 
tions, On the other hand, there are no double entries, either from 
misrecognition or from double radio spots. Though the positions 
unevenness prevent them from characterizing an astrometric ref- 
erence frame, the LQAC positions and magnitudes, are adopted 
as starting point for finding the objects in the dense optical cata- 
logs presented next. Figure 1 presents the LQAC sky distribution 
and sky density. 

2.2. USNOB1.0 

The U.S. Naval Observatory B1.0 catalog (Monet et al., 2003), 
hereafter B 1 .0, is the latest one of the series of all-sky ultra dense 
catalogs issued by the USNO. It brings two epochs positions and 




Fig. 1. Sky distribution (equatorial coordinates, N up, E right) of 
the quasars found in the LQAC catalog. The highest density re- 
gions indicate the SDSS DR5 contribution. The scale represents 
the number of sources in regions of 9 sq° shown on the plots. 
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Fig. 2. Sky distribution (equatorial coordinates, N up, E right) of 
the quasars found in the USNO B 1 .0 catalogs. Again, the highest 
density regions indicate the SDSS DR5 contribution. The scale 
represents the number of sources in regions of 9 sq° shown on 
the plots. Notice that the scale is slightly different from that of 
Fig. 1 relative to the LQAC. 



three colors magnitudes (Johnston B, R and I) for 1,042,618,261 
objects, up to about V=20. It results from scans of 7,435 Schmidt 
plates taken for six sky surveys and calibrated using Tycho-2 
stars (Hog et al., 2000). The nominal astrometric accuracy is 200 
mas, but the individual objects position precision has a median 
value of 120 mas along each direction of measurement. Noticing 
that the quasars forming the LQRF seat on the B1.0 fainter end, 
where astrometric precision and accuracy degrades, the LQAC 
objects were searched for in the B1.0, and a neighborhood of 
stars, i.e, their position, proper motion, and magnitude, was also 
collected around the found objects. Figure 2 presents the sky 
distribution and sky density of the 83,980 quasars recognized in 
the USNO B1.0 
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Fig. 3. Sky distribution (equatorial coordinates, N up, E right) of 
the quasars found in the GSC2.3 catalog. As before, the highest 
density regions mark the SDSS DR5 contribution. The scale rep- 
resents the number of sources in regions of 9 deg 2 shown on the 
plots. Notice that the scale is slightly different from that of Fig. 
2 relative to the B1.0 catalog. 



2.3. GSC2.3 

The Second Generation Guide Star catalog latest version (Lasker 
et al., 2008), hereafter GSC2.3, is an all-sky catalog derived 
from the Digitized Sky Survey (Taff et al., 1990) that the Space 
Telescope Science Institute and the Osservatorio di Torino have 
created from the Palomar and UK Schimdt survey plates. It con- 
tains position, proper motion, and magnitude (Johnston B, R, 
and I) for 945,592,853 objects, and is expected to be complete 
to R=20. The astrometric total error is quoted as smaller than 
300 mas, while the relative astrometric error is better than 200 
mas. In fact the astrometric residuals against external reference 
catalogs only become larger than 150 mas above R=18 at the 
fainter end. The astrometric calibration was done using Tycho-2 
stars. As seen, there are overwhelming similarities between the 
GSC2.3 and the B1.0 raw data, which however were measured 
and treated by distinct methods, enabling a most natural combi- 
nation of the independent results issued from each of them. So, 
as before, the LQAC objects were searched for in the GSC2.3, 
and a neighborhood of stars (position, proper motion, and mag- 
nitude) was collected around the found objects. Figure 3 presents 
the sky distribution and sky density of the 93,943 quasars recog- 
nized in the GSC2.3. 

2.4. SDSS DR5 

The Sloan Digitized Sky Survey (SDSS) data release 5 
(Adelman-McCarthy et al. 2007), hereafter DR5, represents the 
completion of the SDSS-I project. It covers 8,000 sq° , on the 
north Galactic cap (b < 30°) and on three stripes towards the 
south Galactic cap, reaching up to i=20.2, to a total of 215 mil- 
lion unique objects. The astrometry is preferably referred to the 
UCAC catalog, or else to the Tycho-2 catalog. The RMS residu- 
als per coordinate are 45 mas and 75 mas for reductions against 
the UCAC (Zacharias et al., 2000) and the Tycho-2 (Hoeg et al., 
1997), respectively, with additional 30 mas systematic errors in 
both cases (Pier et al., 2003). The internal astrometric precision 
is mainly limited by seeing effects, and it is quoted around 35 



mas in each coordinate. A point to notice is that quasars are sig- 
naled in the DR5, to a total of 74,869 entries with quasi-stellar 
object spectral class. As seen from Fig. 1, the large number of 
DR5 quasars represents an important feature of the LQAC input. 

2.5. UCAC2 

The second USNO CCD Astrograph catalog (Zacharias et al., 
2004), hereafter UCAC2, is the latest release of the ongoing 
UCAC project, designed to observe the entire sky for R mag- 
nitudes between about 7.5 and 16. The observed positional er- 
rors are about 20 mas for the stars in the 10 to 14 magnitude 
range, and about 70 mas at the limiting magnitude of R=16. The 
UCAC2 is a high density, highly accurate, astrometric catalog 
(positions and proper motions) of 48,330,571 stars covering the 
sky from -90 to +40° in declination and going up to +52° in 
some areas. The northern limit is a function of right ascension. 
The astrometry provided in the UCAC2 is on the Hipparcos sys- 
tem. Positions are given at the standard epoch of Julian date 
2000.0. The UCAC2 supersedes the Tycho-2 catalog, having a 
stellar density that is over 20 times higher. It is used here as the 
choice catalog to place the quasar positions collected from the 
catalogs described in the above subsections onto the Hipparcos 
system. Nevertheless, since the UCAC2 coverage does not attain 
the entire sky, two additional stellar frames are complementary 
used and are presented below.. 

2.6. UCACN 

What is hereafter termed UCACN is a cutout of the prelimi- 
nary UCAC positions around the LQAC quasars in the north- 
ern part of the sky that is not reached by the UCAC2 coverage 
(Zacharias, 2007). Although, as discussed further on, the astro- 
metric accuracy of the UCACN positions and its stellar density 
granted no degradation on the results obtained using it, it will 
be dropped in future revisions of this work as soon as of the 
forthcoming release of the all-sky catalog UCAC3. The nominal 
astrometric accuracy of the UCACN stars is similar to that of the 
UCAC2 ones, but no proper motion data were released for these 
stars. Throughout this paper when referring to the UCAC2 and 
UCACN ensemble the designation UCAC is used. 

2.7. 2MASS 

The Two Micron All-Sky Survey point source catalog (Cutri et 
al.,2003), hereafter 2MASS, derives from uniformly scanning 
the entire sky in three near-infrared bands to detect and charac- 
terize point sources brighter than about 1 mJy in each band, with 
signal-to-noise ratio greater than 10, using a pixel size of 2.0", 
each point in the survey having been imaged six times. The de- 
tectors worked to a 3<x limiting sensitivity of J=17.1, H=16.4, 
and K=15.3. The 2MASS contains the position of 470,992,970 
sources, but no proper motions. The astrometry is referred to 
the Tycho-2 catalog and it is accurate to 70-80 mas over the 
magnitude range of 9 < K < 14 mag. Comparative studies have 
shown that the 2MASS positions agree with those from common 
UCAC2 stars to within 10 mas (Zacharias et al., 2003), and that 
the 2MASS is compliant to the ICRF in the range 100-120 mas, 
Taking advantage of its all-sky coverage, very high stellar den- 
sity, and consistent astrometry, the 2MASS catalog is used here 
as an independent vehicle to place the LQAC quasars recognized 
in the large, deep catalogs on the Hipparcos system. 
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2.8. ICRF-Ext2 

The International Celestial Reference Frame 2 nd extension (Ma 
et al.,1998, Fey et al., 2004), hereafter ICRF-Ext2, is the present 
materialization of the International Celestial Reference System 
(ICRS) at radio frequencies. It represents the basic frame with 
respect to which the position of any object in the celestial sphere 
should be measured. In its primary contents, the ICRF consists of 
212 sources, called "defining", whose positions are independent 
of the classical planes (equator, ecliptic) and reference points 
(equinox), but consistent with the previous realizations of the 
Celestial System as the FK5. With the help of a considerable 
amount of VLBI observations, the individual positions of the 
sources were found to be accurate to within roughly 0.25 mas, 
while the stability of the reference axes attains a remarkable 
20 /ias (micro-arcsec) accuracy. In the ICRF initial version 608 
sources complemented the defining 212 sources. The entrance of 
the ICRF-Extl and of the ICRF-Ext2 brought 109 new sources in 
the catalog, thus leading to a total number of 717 radio-sources. 
Notice that in the ICRF-Ext.2 catalog, a very small sample of 
objects are not quasars: 10 of them are cathegorized as AGN 
(Active Galactic Nuclei) and 10 of them are cathegorized as BL 
LAC (BL Lacertae). Although the LQAC a priori exclusively 
considers quasars for its compilation, these particular objects 
were retained thanks to their astrometric accuracy. The ICRF- 
Ext2 is used here to match the VLBI positions directly to the de- 
rived optical positions. Such match serves to two distinct steps: 

(a) to redress the global frame orientation towards the ICRF, and 

(b) to redress local departures of the derived quasar catalog from 
the ICRF. 

2.9. VCS6 

The Very Long Baseline Array (VLBA) Calibrator Survey 
(Petrov et al., 2006), hereafter VCS6, consists of a catalog con- 
taining milli-arcsec accurate positions of 3,910 extragalacic ra- 
dio sources, mainly quasars. These positions have been de- 
rived from astrometric analysis of dual-frequency 2.3 and 8.4 
GHz VLBA observations, on the framework of the International 
Celestial Reference System (ICRS) as realized by the ICRF. 
Thus, to the level of precision and accuracy of the optical posi- 
tions determined here, the VCS6 sources represent the ICRF. In 
order to strictly validate such representation, 524 VCS6 sources 
with formal inflated errors larger than 10 mas were trimmed off. 
The VCS6 sources not belonging to the ICRF-Ext2 were added 
to it, to obtain an enlarged radio frame, which is used for the two 
steps mentioned in the ICRF-Ext2 subsection, that is, the global 
and local redressing of the optical quasar positions towards the 
ICRF. 

2.10. VLAC 

The Very Large Array (VLA) Calibrator list (Claussen, 2006), 
hereafter VLAC, consists of 5,523 radio flux, structure code, 
position, and associated error entries for 1,860 extragalactic 
sources. It is foremostly indicated as a reference tool for ob- 
servations at the VLA, but the radio positions are accurate and 
compliant with the ICRF. Then, likewise in the previous sub- 
section, the VLAC positions were added to the ICRF-Ext2 core 
data set, after trimming off 297 of them with errors larger than 10 
mas. The combination of the ICRF-Ext2, VCS6, and VLAC will 
be here termed as the enlarged radio frame (ERF) and it repre- 
sents the ICRF to within an uncertainty of 10 mas. The order of 
priority for the radio lists to compound the enlarged radio frame 




Fig. 4. Sky distribution (equatorial coordinates, N up, E right) of 
the enlarged radio frame quasars. The filled symbols mark the 
sources for which an optical counterpart is found in at least one 
of the deep, dense source catalogs (B1.0, GSC2.3, DR5). The 
enlarged radio frame is formed, in order of priority of choice, by 
the ICRF-Ext2 (718 sources), the VCS6 (2,684 sources), and the 
VLAC (59 sources). 

was, ICRF-Ext2 (718 sources), VCS6 (2,684 sources), and fi- 
nally VLAC (59 sources). The 3,461 sources now represent the 
ICRF on a better distributed all-sky basis (Fig. 4), and those for 
which the optical counterpart position is determined will actu- 
ally define the best adherence of the LQRF to the ICRF. 

2.11. Data flow 

The various origins of the sources compiled in the LQAC im- 
poses first of all to refer their positions to a single frame in the 
ICRS system. After that, the quasar frame so formed itself is 
aligned with the ICRF. Finally, the position of each source is 
calculated by weighted average. This data flow can be sketched 
as follows. 

1. The LQAC entries are all admitted. 

2. The LQAC quasars are searched for, by a main criterion of 
position coincidence, in the B1.0, GSC2.3, and SDSS DR5 
high density catalogs. Each individual match is kept and 
treated separately. 

3. Around each of the individual quasars thus matched a stel- 
lar neighborhood is extracted from the pertinent high density 
catalogs. 

4. The stellar neighborhood contents from the high density 
catalogs are searched for in the stellar catalogs UCAC2, 
UCACN, and 2MASS, as before by a main criterion of posi- 
tion coincidence. 

5. Within each neighborhood, the stars common to the high 
density catalog and to the stellar catalog have their positions 
projected onto the tangential plane centered at the quasar. A 
local transformation function from the high density catalog 
to the stellar catalog is calculated, mimicking the complete 
plate solution polynomials. 

6. At this point, for each quasar, typically, a number of posi- 
tions are determined (e.g., from more than one high density 
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catalog, and local solutions from the two stellar catalogs). 
Each of these families of positions will still be handled sep- 
arately. 

7. Within each family, the quasars for which there is a position 
in the enlarged radio frame (that is, the combination of the 
ICRF, VCS6, and VLAC catalogs) are found, since they are 
so flagged in the LQAC, what as before follows mainly a 
criterion of positional coincidence. 

8. The subsets so defined are used to calculate the global ro- 
tation and zero point corrections to obtain coincidence with 
the ICRF origin and axes. The global rotation and zero point 
corrections are next applied to the positions of all quasars of 
each given family. 

9. Likewise the previous step, the optical and radio positions in 
subsets are combined using orthogonal functions, in right as- 
cension, declination, and magnitude, in order to correct sys- 
tematic local departures from the ICRF axes directions. And 
likewise the significant orthogonal functions are applied to 
the positions of all quasars of each given family. 

10. Again the optical and radio positions in subsets are combined 
in small regions around each quasar to derive its correction 
for localized inhomogeneities. 

11. For each position of each source in a given family, a total 
error is assigned that combines the formal error from the high 
density catalog where it appears, the error from the particular 
local correction, the error from the particular global rotation 
and bias, and the error from both the significant orthogonal 
function adjustments and the inhomogeneity corrections. 

12. From the squared inverse of the total error, a weight is as- 
signed to each determined position of a given quasar. The 
final position of each quasar is given as the weighted aver- 
age of the determined positions. The error in the position is 
assigned according to an individualized function of the opti- 
cal to radio disagreement. 

The comparisons of the final positions against the existing 
corresponding positions in the enlarged frame is used to obtain 
the global properties of the LQRF. 

3. The local astrometric solution 

In this section the points 1 to 6 of Sect. 2.11 (data flow) are 
derived. 

Starting from the LQAC input list of quasars, the first step to 
get their astrometric positions, in order to form a whole consis- 
tent reference frame, is to have them matched to the position of 
objects in the two large, deep all sky catalogs, the USNO B1.0 
and the GSC2.3, as well as in the SDSS DR5. The basic match- 
ing criteria is of position agreement within 1 arcsec. This is equal 
to the poorest astrometric accuracy in the LQAC contributing 
surveys. It is also a safe threshold given the typical seeing and as- 
trographic plate resolution limits. Additionally a magnitude limit 
was imposed preventing the selection of objects brighter than 
R=7, nevertheless this threshold was not reached by any source. 
No proper motion limit was enforced because both the USNO 
B1.0 and the GSC2.3 warn against possible offsets of the proper 
motions zero points. A post check was made by determining lo- 
cal zero points of proper motion, given by the average motion. 
In comparison to the local zeros of proper motion, just 0.7% of 
the GSC2.3 selected objects show significant motion, while for 
the B 1 .0 selected objects just 0. 1 % of them show significant mo- 
tion. To verify what would be the possibility of a chance hit, the 
input coordinates were varied by random values between 1 and 5 
arcmin. In this way the same region of the true objects was being 
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Table 1. Initial optical to radio relationships 3 . 



Tracer Bl.OtoRadPos GSC2.3 to RadPos DR5 to RadPos 



N 


1962 


2028 


300 


AacosS 


-35.2 ±4.7 


+32.9 ± 5.5 


-4.5 + 3.2 


A6 


+ 123.5+4.9 


+ 39.5 ± 4.7 


+ 13.9 + 3.2 


Al 


+4.7 ± 6.2 


-6.6 ± 6.7 


-3.5 ± 4.7 


A2 


+5.9 ±5.8 


+ 13.1+6.3 


-19.4 + 5.1 


A3 


+42.1 ± 5.7 


-29.0 + 6.1 


+6.6 + 4.1 


A4N 


+ 123.0 + 5.3 


20.9 ± 5.7 


1.6 + 4.5 


A4S 


+ 128.0 ± 11.9 


134.1 ± 12.5 





a For the catalogs USNO B 1 .0, GSC2.3, and SDSS DR5 relatively to 
the collected sample of sources with precise radio interferometric 
positions (Extended Radio Frame - RadPos, either from ICRF-Ext2, 
VCS6 or VLAC), are presented: the average offsets in right ascen- 
sion and declination, the orientation angles (Al, A2, A3) relative to 
the standard thriedron of equatorial coordinates, and the equatorial 
biases (A4N and A4S, north and south of -20°). All values are in 
milli-arcsec (mas) and the number of sources is given in the first 
line. 



swept. Keeping the search radius to 1 arcsec, in only 0.2% of the 
cases was verified a hit on the false coordinates. Moreover, in 
only 0.1% of cases was a second object found significantly close 
to the adopted match. These satisfying proportions change very 
little when the search radius is enlarged to 2 arcsec. While this is 
an added proof of the overwhelming correctness of the matches, 
it at the same time indicates that there is little advantage in adopt- 
ing a larger radius, and certainly no gain to the confidence level 
attached to the cross-identifications. 

In all, 83,980 quasars are recognized in the B1.0, and 93,943 
quasars are recognized in the GSC2.3. The sky distribution fol- 
lows closely that shown in Fig. 1. From the figure it is clear 
the large contribution of the SDSS DR5. For the later, 74,825 
quasars are retrieved using the SDSS query server. The positions 
from the three catalogs were added independently to the input 
list. The final input data contains 100,165 quasars. Just small 
fractions of the total appear in only one of the catalogs (1,405 
only in USNO B1.0, 4,425 only in GSC2.3, and 3,276 only in 
SDSS DR5), while 61,532 quasars appear in all three catalogs. 

These three catalogs that bring the input positions for the 
recognized quasars are originally placed on the ICRF J2000 
reference frame by different astrometric pipelines based on the 
Tycho-2 catalog. Their global orientation and equatorial bias rel- 
ative to the ICRF J2000 frame are shown in Table 1 . Local devi- 
ations are also apparent in Figs. 5-7. 

In order to homogenize all the original different positions, 
a local astrometric reduction was performed based on the cata- 
logs UCAC and 2MASS. The local correction follows the stan- 
dard tangent plane astrometric reduction, taking the quasar orig- 
inal position as the central point of the rigorous gnomic projec- 
tion (Assafin, Vieira Martins & Andrei, 1997). Through it the 
catalogs coordinates are adjusted to the reference frames repre- 
sented by the stellar catalogs. Three assumptions are implicit to 
the method. That the stellar catalogs deliver an improved, and 
similar, representation of the ICRF, either relatively to the two 
Schmidt plate surveys, as well as relatively to the photometric 
one. That in a very restricted neighborhood the deviation be- 
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Fig. 5. USNO B 1 .0 right ascension and declination offsets distributions relative to the equatorial coordinates axes and the R magni- 
tude. The offsets are given in the sense of catalog minus the radio interferometric positions from the ICRF-Ext2, or the VCS6, or 
the VLAC. 



tween the deep surveys and the stellar catalogs, whichever one it 
is originated from, can be modeled by simple relationships. And 
that the removal of such deviations can significantly improve the 
accuracy of the quasar positions found in the large surveys. Prior 
investigations support these assumptions (da Silva Neto et al., 
2005). 

3.1. The polynomial degree of the solution 

On the functional side, the initial question to be tackled is on 
the polynomial degree of the local astrometric correction repre- 
sentation. The UCAC is the least dense of the involved catalogs, 
and it presents a typical density of 10 stars per 5x5 arcmin. 
Translated to the Schmidt plates scale this means that a first de- 
gree solution with 10 reference stars can be applied over a phys- 
ical area as small as 25mm 2 . According to the adopted a priori 
assumptions, and given the sky density of the stellar catalogs, 
the simplest models seem to be well suited. To verify such pro- 
cedure, a trial local reduction was performed for 29 quasars ran- 
domly distributed across the celestial sphere, at every one hour 
of right ascension, off the galaxy plane. Around each of them, a 
field of UCAC2 stars was selected in boxes of side 10 arcmin, 
which contained at least 30 stars. For these 29 fields of quasars, 



complete polynomial solutions were tested, with orders varying 
from to 3 rd , and a straight averages solution. The quasar po- 
sitions were taken from the USNO B 1 .0 and all positions were 
transported to the USNO B 1 .0 plate epoch, using the proper mo- 
tions of either catalog. Table 2 summarizes the mean results. 

In the first half of Table 2 the average local corrections and 
the corresponding standard deviation are shown. Though the val- 
ues only very loosely represent the systematic offsets between 
the USNO B 1 .0 and the UCAC2, what is important to note is the 
similarity between all the rows. In the second half of the table, 
the average internal errors of each solution and their dispersion 
are shown. Again all rows behave quite alike. Therefore, to keep 
the smallest the region where the local solution is derived, and 
to maximize the ratio of the number of equations to the num- 
ber of parameters to find, the complete first degree polynomial 
was adopted. This solution also produces the smallest set of stan- 
dard deviations. The zero degree solutions fared just as well as 
the others, so there would be no particular advantage on taking 
an even smaller region, while the first degree polynomial allows 
for some modeling, which is useful face to the large numbers of 
sources that the local solution is going to be applied to. Two for- 
mulations of the first degree polynomial are concurrently used, 
the complete one (six independent parameters) - hereafter the six 
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Fig. 6. GSC2.3 right ascension and declination offsets distributions relative to the equatorial coordinates axes and the R magnitude. 
The offsets are given in the sense of catalog minus the radio interferometric positions from the ICRF-Ext2, or the VCS6, or the 
VLAC. 



parameter model, and the so-called four constants one (in which 
the X and Y axes are constrained to have the same scale) - here- 
after the four parameter model. Furthermore, notice that with the 
six parameter model the error in the central, pseudo-tangential 
point, which is always the position of the quasar, was also found 
to be the smallest. 



The 2MASS Point Source catalog positions are compliant 
with the ICRS via the Tycho-2 catalog and expected to be accu- 
rate to 70-80 mas over the magnitude range of 9 < Ks < 14 mag. 
At the limit Ks = 16, the accuracy is not better than 200 mas. As 
the 2MASS does not provide proper motions, these values are 
good for the mean epoch of observation, 1999. 



3.2. The stellar catalogs and the size of the local segion 

As discussed above, the UCAC2 offers the desiderata required 
for the local astrometric solution: high stellar density, accurate 
global representation of the HCRF and hence of the ICRF, range 
of magnitudes up to 16, and small zonal and local frame biases. 
The stellar density supports using a region as small as 10 arcmin 
across, which is always within the overlap zone of the Schmidt 
plate surveys contributing to the B1.0 and the GSC2.3. 

However, the UCAC2 has incomplete coverage in the north- 
ern celestial hemisphere. Two remedies are used to face such 
hindrance, the use of preliminary positions of northern UCAC 
stars (here called UCACN), and the use of the 2MASS catalog. 
In the former case, the UCACN positions are by construction on 
the same system of the UCAC2, and at the same precision. 



The pitfalls due to the magnitude gap between the 2MASS 
objects and the brighter Tycho-2 reference stars degrade the fi- 
nal astrometric accuracy. The internal precision of the 2MASS 
individual positions is better, at the level of 40-50 mas in the 
9 < Ks < 14 mag range. The comparison between 5.2 million 
stars common to the UCAC2 shows standard deviations of 70- 
80 mas in both coordinates, for the 9 < Ks < 14 mag range. For 
the common stars around the quasars extracted to define the lo- 
cal solution, the difference average and standard deviation are 
Aacosd = +3 mas, cr AncosS = 99 mas, and A6 = +11 mas, cr M 
= 95 mas, all in the sense 2MASS. minus UCAC2. Furthermore, 
from the limited sample of 700 quasars from the enlarged radio 
frame found in the 2MASS catalog, the analogous statistics are 
Aacosd = +0 mas, cr &acosS = 144 mas, and A6 = -2 mas, o~as = 
138 mas. 
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Fig. 7. SDSS DR5 right ascension and declination offsets distributions relative to the equatorial coordinates axes and the R magni- 
tude. The offsets are given in the sense of catalog minus the radio interferometric positions from the ICRF-Ext2, or the VCS6, or 
the VLAC. 



On the positive side, for the same region, the 2MASS catalog 
presents about 10 times the number of reference stars that does 
the UCAC catalog. In the small regions used for the local solu- 
tion, it represents a coherent local frame, from which the zonal 
bias can be removed after later comparisons against the enlarged 
radio frame sources. It is well aligned with the UCAC and with 
the ICRF, and represents an all-sky frame. The local solution 
was computed independently using UCAC and 2MASS refer- 
ence stars. Each of the results contribute, again independently, 
to the weighted average used to compute the final LQRF quasar 
coordinates. 



3.3. The analytical expression for the local correction 
astrometric solution 

On establishing the analytical form of the local astrometric so- 
lution, a final point to be considered concerns the usefulness of 
a magnitude dependent term, a feature often present in Schmidt 
plates based surveys (da Silva Neto et al., 2000). In the present 
case, the R magnitude represents the obvious choice, mainly be- 
cause it is common to all the involved catalogs, and closer to 



the central wavelengths of the 2MASS bands. At the same time, 
due to the quasars redshift and to the reddening for the stars, 
coupled with CCD top sensitivity and wider plate effective pass- 
band, the signal-to-noise ratio is usually higher towards the red. 
An inspection on the magnitude dependency of the catalog mi- 
nus enlarged radio frame residuals reveals a small linear trend, 
larger for the A6 offsets than for the AacosS offsets. The same 
trend is verified for both B1.0 and GSC2,3 increasing with mag- 
nitude for declination, and decreasing with magnitude for right 
ascension. The amplitudes get only to 10 mas in the significant 
cases. In what concerns the DR5, the magnitude dependence is 
of much smaller amplitude and opposite tendency. The possibil- 
ity of a magnitude equation is not ruled out in the presentation of 
either of the two deep catalogs. However, the causes and ensuing 
description would be out of the scope of this work. Here, an ad- 
ditional obstacle would be to bridge the magnitude gap between 
the reference stars and the much dimmer quasars. All accounted, 
the magnitude investigation is left to be investigated by the or- 
thogonal functions in direct comparison with the enlarged radio 
frame positions. 

Therefore, the enforced local solutions solution are: 
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Table 2. Comparison between the models for the local astromet- 
ric solution 3 . 



Solution 


Aacosd 


C Aacosd 


AS 


o-as 


yd dggj-gg 


+42 


126 


-25 


140 


2 nd degree 


+43 


127 


-23 


136 


1" degree 


+61 


110 


-34 


120 


degree 


+42 


126 


-34 


119 


offset 


+56 


118 


-45 


130 



Solution 






S r 


0--LY 


yd (Jggj-gg 


137 


63 


127 


64 


2 nd degree 


134 


66 


125 


64 


I s ' degree 


135 


68 


120 


65 


degree 


133 


68 


117 


64 



In the top half of the table, the average corrections and their stan- 
dard deviations are presented. In the bottom half, the average in- 
ternal error in the solutions and their dispersion are presented. All 
values are in milli-arcsecs. The difference between the "0 degree" 
and the "offset" solutions is just that in the former residuals larger 
than 2.5<t are removed and the solution is recalculated. 



%SCat ~ XQCat - a X Q Cat + b YqCci + c (1) 
f] 'SCat ~ YQCat = d Yqcm + e X Q c at + f (2) 



In Eqs (1) and (2), "SCat" is either the UCAC2 (or even its 
preliminary northernmost part) or the 2MASS, and "QCat" des- 
ignates the B1.0, the GSC2.3 or the DR5. The relations above 
are fit to the smallest set of common "SCat" and "QCat" stellar 
objects in a neighborhood larger than 10 arcmin but smaller than 
a box of size of 30 arcmin, to a minimum of 6 stars. The solu- 
tions are calculated independently for right ascension and dec- 
lination. Whenever that minimum number of stars is touched, 
either at the start or after elimination of stars, the four param- 
eter model is used - in that case, obviously combining the two 
equatorial coordinates solutions. In fact, as presented in the next 
section, the four parameter model was run for all quasars and it 
delivers solutions in excellent coincidence with the six parame- 
ter solutions. The analytical form of the four parameter model is 
presented in Eqs. (3) and (4), where "SCat" and "QCAT" keep 
their afore-defined meaning. 

%SCat ~ Xqcm = gXQCat + h YQCat + i (3) 
rjSCat ~ Y Q Cat = k Y QC at - gX QCa t + j (4) 



3.4. Applying the local astrometric solution 

As discussed in the previous sections, the initial task of produc- 
ing the LQRF consists of carrying out a local astrometric reduc- 



tion using the most precise and accurate stellar catalogs repre- 
senting the HCRF, through the Tycho-2 reference frame, in or- 
der to locally redress the quasars' position input from the high 
density, deep catalogs.. 

There, however, the question of epochs comes up. The stel- 
lar catalogs have a mean epoch very close to J2000, being 
J1998.79 (+ 2.78 y) for the UCAC2, and J2003.34 (± 0.24 y) 
for the UCACN, the UCAC2 preliminary northernmost part, and 
J1998.87 (+ 0.08 y) for the 2MASS. Also, the DR5 mean epoch 
is J2002.96 (± 1.45 y), raising no questions. The dilemma on 
whether the local reduction ought to be performed in the epoch 
of the quasar input position or directly in J2000 is actually only 
meaningful for the cases when the B1.0 and the GSC2.3 are lo- 
cally reduced by the UCAC2 . For them it would not be unrea- 
sonable to speculate that a reduction made at the quasar epoch 
would diminish the scatter in the stellar positions, on the suppo- 
sition that the stellar images themselves were at the quasar epoch 
(i.e., the epoch of the plate where the quasar is found) because of 
the smallness of the used field. The B1.0 mean epoch is J 1977. 12 
(± 5.49 y). On the other hand, it is verified that the B1.0 posi- 
tions are given as an average of plates reductions, in such a way 
to jeopardize the hypothesis of common epochs, given the range 
of magnitudes and to some extent the two different types of ob- 
jects. For the GSC2.3, in contrast, the positions always refer to 
one single highest quality plate, whenever possible the red one. 
But, again, the GSC2.3 mean epoch is J1992.5 (+ 3.72 y) and 
therefore there is a balance between reducing modestly the stel- 
lar positions scatter and the errors brought in by faulty or biased 
proper motions. 

In order to actually verify which way would be preferred, the 
local corrections for the sources belonging to the enlarged ra- 
dio frame were evaluated both at J2000 and the quasar position 
epoch. A second trial was made by considering only the sources 
from the enlarged radio frame belonging to the three input cat- 
alogs and that could be reduced by all the stellar frames. This 
subsample contains 240 sources but enables a direct comparison 
between all the possible ways of treatment. The outcomes show 
no net advantage in working at the quasar epoch. For the B1.0 
catalog, the number of sources for which the local reduction con- 
verges drops; there is a marginal gain on the number of stars 
entering in the reduction but the number of actually used stars 
remains the same; there is no diminishing on the standard devia- 
tions from the solution corrections, or from the optical minus ra- 
dio offsets; and the removal of the equatorial bias is less efficient 
(by 98 mas). The external check actually shows a small increase 
of the average detachment to the representation of the ICRF ori- 
gin. For the GSC2.3 there is no important loss in working at 
the quasar epoch, but no gain either. The number of sources for 
which the local reduction converges remains the same, and so 
does the number of entry and used stars. The standard devia- 
tion in the solution does not change, and for the optical minus 
radio positions the offset scatter improves only marginally (by 
4.5 mas). The equatorial bias is less efficiently removed (by 40 
mas). Exactly the same results are obtained when inspecting the 
subsample of common sources. In view of these results it was 
decided to work only at J2000 for the local corrections. Thus 
all pairs of input quasar catalogs and stellar catalogs are treated 
identically. 

A second point concerning the issue of epochs has already 
been touched in the previous discussion. It concerns whether the 
local reductions made with the UCACN, that is for the north- 
ernmost sources, and with the 2MASS would be significantly 
worse than those made with the UCAC2. As seen, both of them 
have a mean epoch close to J2000, where the local corrections 
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are performed. For the 2MASS catalog, which is not originally 
presented as an astrometric frame, the number of equations of 
condition is substantially increased for a same patch of the skies 
as defined for the UCAC2. To test the performance of the differ- 
ent catalogs, we use the DR5, which was shown in Sect. 3 to have 
intrinsically smaller scatter relative to the ICRF. As before, the 
enlarged radio frame and its subset of all common sources are 
employed, and they provide responses much alike. The UCAC2 
based reductions produce the best results as expected. The right 
ascension and declination standard deviation in the offsets to the 
radio frame is at 47.5 mas. We remark that this is already a gain 
relative to the DR5 catalog offsets that are at 59 mas, as produced 
from Table 1 values. The standard deviation in the solution itself 
is at 33 mas. The UCACN based reductions perform much alike, 
being actually slightly better on declination and 30 mas worse on 
right ascension. The optical minus radio standard deviation is 14 
mas worse than that of the UCAC2. The 2MASS based reduc- 
tions fare somewhat worse than the UCACN ones, being larger 
by 18 mas in the solution standard deviation, and worse by 6.5 
mas in the optical minus radio standard deviations. In this case, 
it is worth to point out that the average optical to radio position 
offsets are much similar for the UCAC2 and the 2MASS (dif- 
fering by just 5 mas). It is also to be stressed that, eventhough 
not bringing a direct improvement in accuracy as the UCAC2 
based ones do, the UCACN and the 2MASS local reductions en- 
able to place the DR5 positions on the same frame with the B 1 .0 
and GSC2.3 reduced positions, without impairing their intrinsic 
precision. 

Table 3 summons up all the applied local corrections. 

As can be seen from Table 3, the different quasar input cata- 
logs (codes B, G, or S) require different local corrections, while 
the three stellar reference catalogs (codes U, N, or T) deliver 
such corrections at similar accuracy. The local correction out- 
puts (codes C) agree when calculated either by the six parameter 
or by the four parameter models (codes 6 or 4). The local cor- 
rections are generally larger than the remaining optical minus 
radio offsets (codes A), which supports their usefulness. The op- 
tical minus radio mean offset in all cases is improved relatively 
to the one obtained using original quasar catalogs positions. As 
a consequence, the mean offsets become similarly small. 

At this point, the three quasar input catalogs and the three 
stellar reference frames deliver concurrent, independent sets of 
positions. And whenever the 1 " degree polynomial solution fails 
to converge, the four parameter solution results can be adopted 
equivalently. There are thus up to twelve possible local solu- 
tions calculated for each source. For 61% of the sources more 
than one solution was obtained, for 844 sources all twelve so- 
lutions went through, and for 12 sources no solution converged. 
Figures 8 to 1 1 display the local corrections for the most critical 
cases, namely the B1.0 and the GSC2.3, as given by the UCAC 
and 2MASS frames. Comparing with the Figs. 5-6, it is already 
clear that the local corrections bring the quasar's positions onto 
the ICRF. The comparison between these figures and the initial 
plots (Figs. 5-7) show that the local astrometric solution signifi- 
cantly redress the input positions from the large catalogs. Notice 
that the outcomes from the UCAC and 2MASS stellar reference 
frames bear much in common. Likewise, the six parameter and 
the four parameter models derive coincident corrections (thus, 
only those from the six parameter model are shown in the fig- 
ures). It is also interesting to note that there is no discontinuity 
between the zones where either the UCAC2 or the UCACN stel- 
lar frames were used. 



4. The global orientation towards the ICRS 

In this section the points 7 and 8 of the Sect. 2.11 (data flow) are 
derived. 

Once the local correction is obtained, the next move is to 
tie the preliminary frame just derived to the ICRF. As discussed 
in Sect. 1, the ICRF alone is too sparse, and in particular the 
number of its bright optical counterparts, to suffice for a detailed 
tying. Taking advantage of the fact that the long base radio inter- 
ferometric observations always combine data for several sources 
in different sessions, a radio position net considerably larger than 
the ICRF was formed by adding the VLBA calibrators and the 
VLA calibrators. In the case of the VLA, only sources with the 
highest astrometric precision are included. The accurate, ICRF 
representing, radio interferometric positions of the enlarged ra- 
dio sample sources (ERF) are then compared with their opti- 
cal counterpart from the stellar based preliminary frame derived 
from the local corrections. From the ICRF 718 sources are col- 
lected, while from the VSC6 and the VLAC, 2,684 sources and 
59 sources are collected respectively. 

The enlarged radio frame so gathered thus contains 3,461 
sources. They are adequately represented in all of the three con- 
tributing catalogs as detailed in Table 3 (column Nrad). The sky 
distribution was shown in Fig. 4. From the ERF sources, 2,263 
have an optical counterpart in at least one of the input catalogs. 
The average minimum distance between the ERF sources is 187 
arcmin, with a mode of 70 arcmin. 

The ERF radio positions are then compared with their optical 
counterpart from the stellar based preliminary frame derived by 
the local corrections. Initially the comparison entails the deriva- 
tion of global orientation corrections about the equatorial coordi- 
nates standard triad, plus a correction for a bias relatively to the 
equator. Following (Arias, Feissel & Lestrade, 1988), the canon- 
ical equations relating the right ascension and declination offsets 
to the rotations (Al, A2, A3) about the equatorial axes (namely, 
X on the equator pointing to the conventional origin, Y on the 
equator perpendicular to X, and Z pointing to the conventional 
pole) are given by 

Aa cos8 = A 1 sinS cosa + A2 sinS sina - A3 cosd (5) 
A(5 = -Al sina + A2 cosa + A4 (6) 



Equations (5) and (6) are solved independently as a mea- 
sure to assess the robustness of the corrections. In this case, the 
direction cosines Al and A2 agree when calculated either from 
Aacos(5 or A6 within 1 cr in most cases. Figure 12 shows the good 
agreement on Al and A2 as calculated from the right ascension 
residuals, from the declination residuals, or from the combined 
residuals. The equatorial bias (which is discussed in detail in the 
next section) can also be independently derived as simply A6 = 
A4. Also for A4, either the straight average or the adopted com- 
bined solution agree well. 

For the sake of gathering the largest number of equations 
of condition, the final values are derived by a combined matrix 
including the three rotation corrections, plus the representation 
of the equatorial bias. 

4. 1 . The equatorial bias 

Former investigations (Assafin et al., 2001) pointed out marked 
differences in the equatorial bias concerning the USNO A2.0 cat- 
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Comparative results given by all types of local correction designed. The solutions nomenclature is formed by one letter designating the quasar 
position input catalog, followed by one letter designating the stellar reference frame, and by the codes 6 and 4 informing whether the algorithm 
is of the six or four parameter model. Nq is the number of quasars for which a given local correction was possible and succeeded. Ns is the 
average number of catalog stars contributing to the solution. CacosS and CS columns bring the mean corrections on right ascension and 
declination, and the standard deviations (in brackets). The final two columns show the average optical to radio offset, preceded by the number 
of sources entering the average (Nrad). All angular values are in milli-arcsec. See the text for detailed explanations on the solutions and 
samples. 



cally corrected by the UCAC2 and the 2MASS frames, adopting 
a single A4 unknown, two A4 unknowns dividing the AS offsets 
at -20° of declination, and eighteen A4 unknowns to fit the AS 
offsets within declination strips 10° wide. The plots in Fig. 13 
show features that are more complex than could be described 
by one A4 standard offset. Nonetheless, there is a coarse posi- 
tive plateau southwards of declination -20°, and values averag- 
ing near zero northwards. The standard deviation in the residuals 
from the 18 strips solutions is indeed no smaller than that of the 
-20° boundary, 2 A4 terms solution. This testifies of an inter- 
play that can be more adequately handled in the next step where 
harmonic terms are introduced. On the other hand, the solution 
contemplating a single A4 term fares the worst by a few mas, 
and is certain to provide the least appropriate description. 

The verified A4 declination dependency, prompted to assess 
whether a similar right ascension dependency is present. In this 
case, no such dependency exists. 

4.2. Rotation angles and equatorial bias 

The final values of the Al, A2, and A3 direction cosines for the 
global orientation towards the ICRF, as well as the equatorial 
bias, were found by adjusting the optical minus radio residual 
between the positions evaluated from the local astrometric solu- 
tion and those given by the enlarged radio frame. The values are 
presented in Table 4. 

For the B1.0 and GSC2.3 northernmost sources, for which 
the local correction was made using the UCACN, separate val- 
ues were derived for the rotation angle and the (northern only) 
equatorial bias. In what concerns the DR5 northernmost sources, 
though the number of UCACN treated sources is small, the pre- 
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Fig. 12. Coincidence of the Al and A2 direction cosines as ob- 
tained either from the AacosS residuals, from the AS residuals, 
or from both. The value adopted was from the solution including 
both the right ascension and declination residuals. 



alog. The point of discontinuity seats around -20° of declination, 
which is the boundary of the north and south plate surveys that 
contributed to that catalog. Although such marked discontinuity 
should not befall neither the B 1.0 nor the GSC23 catalogs, they 
mostly share the same constituent surveys with the USNO A2.0. 

In order to analyze whether a discontinuity remains in the 
data, Eq. (6) was solved for the catalogs B 1.0 and GSC23, as lo- 
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Fig. 8. Equatorial coordinates distribution of the local astrometric corrections derived by the six constants model using the UCAC 
stellar reference frame onto the USNO B1.0 quasar positions. The clear points correspond to the quasars for which there are ERF 
radio positions. 



cisions are high and so a complete solution including one equa- 
torial bias term was sufficient. At large, for the DR5 the Al and 
A2 values, although significant and representing true corrections 
to be applied to the quasar equatorial coordinates, must be inter- 
preted with caution in terms of what they would represent for the 
DR5 frame orientation, because of the concentration of the DR5 
sources in the direction of the right ascension of 12h. 



5. The removal of zonal bias 

In this section the points 9 and 10 of the Sect. 2.11 (data flow) 
are derived. 



up the axes from waves of corrugation either trapped in the stel- 
lar catalogs themselves, or that were present in the initial deep 
catalogs and could not be completely removed by the local cor- 
rections. 

The problem, the quantity of assessing points, the sky cover- 
age, the dimensions of dependencies, and the relationship from 
the less precise, corrugated frame to the more precise, correct- 
ing one, are all the same to the case treated by Schwan (1988) 
when deriving zonal corrections to the FK5 construction. Thus, 
the same treatment by spherical harmonics (Brosche, 1966) is 
used here. The right ascension and declination corrections are 
expressed by the linear combination of the terms of a series, in 
which the cpnml - th element is given by 



5.1. Harmonic functions 

At this point, the optical frame was homogenized by the local 
corrections, and was globally aligned with the ICRF. The final 
point is to sweep its axes to check for any remaining zonal modu- 
lations. This was done for the right ascension, declination, and R 
magnitude dimensions. In a sense this is equivalent to straighten 



Ac = C cpnml Hp(R') Ln(sinS) Fml(a) 



(7) 



where Ac is one or the other equatorial coordinates offsets, 
since independent spherical harmonic series are derived for each 
of them; Ccpnml is the cprml - th coefficient to be derived, for 
the c equatorial coordinate series; Hp(R') is the Hermite poly- 
nomial of degree p, function of the R magnitude normalized by 
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Fig. 9. Equatorial coordinates distribution of the local astrometric corrections derived by the six constants model using the 2MASS 
stellar reference frame onto the USNO B1.0 quasar positions. The clear points correspond to the quasars for which there are ERF 
radio positions. 



the mean and standard deviation for the quasars; Ln(sinS) is the 
Legendre polynomial of degree n, function of sinS; and Fml(a) 
is the Fourier term of order m on a and sign I. 

The offsets are produced by comparing the equatorial coor- 
dinates from the previous step (global rotation) with the VLBI 
positions from the ERF. That is, 2.263 pairs are available. Again, 
this is similar to the quantity available for the Schwan's (1988) 
application. Therefore terms up to very higher order can be de- 
rived. However, the analysis from the earlier steps advises to 
some caution in defining the actually significant order up which 
to halt. The distribution in right ascension of the offsets from the 
original deep catalogs hints for a 2n term and smaller ripples (cf. 
Figs. 5-7), while the analysis of global rotation coefficients has 
shown no a dependence. In contrast, the analysis of the equato- 
rial bias displays a 6 dependence, coherent in bins of 10° (Fig. 
13). Magnitude dependences are believed to be possible within 
the large quasar catalogs. They could to some extent be intro- 
duced by the local solution itself, because of the magnitude dif- 
ference between stars and quasars. The analysis made in Sect. 
3.3 points out to such dependences but of feeble amplitude. All 
accounted for, and in order to have representative quantities of 



ERF objects with which to compute the coefficients of each term, 
it was adopted to extend to 9 th order the Lagrange {sinS) polyno- 
mials and the Fourier (a) functions. The Hermite (R magnitude) 
polynomials were extended to 2 nd order, but just to the 5 th order 
of the other two terms. Given the poor sky distribution, no dec- 
lination terms were used for catalog combinations including the 
UCACN frame, nor right ascension terms for catalog combina- 
tions including the DR5 quasar list. All the above choices were 
made to retain only truly significant terms derived from the ERF 
comparison, such that they could be confidently applied to the 
full set of quasars. 

Since by definition the spherical harmonics are orthogonally 
independent, their statistical significance can be tested individu- 
ally. Each term was initially tested separately. Only those statis- 
tically significant at 3<x with coefficients not smaller than 1 mas 
were retained. Once more, this taxes heavier on the high order 
polynomials, Next, starting from the lower orders and up, the 
corrections brought by each term are removed from the equa- 
torial offsets, and the following term is individually tested. The 
same criterion of 3<x significance and coefficient not smaller than 
1 mas is applied to retain the terms. Actually, at this round the 
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Fig. 10. Equatorial coordinates distribution of the local astrometric corrections derived by the six parameter model using the UCAC 
stellar reference frame onto the GSC2.3 quasar positions. The clear points correspond to the quasars for which there are ERF radio 
positions. 



uppest order Hermite magnitude polynomials, when convolved 
with the higher Legendre (declination) and Fourier (right ascen- 
sion) orders, are retained only if the coefficient is not smaller 
than 5 mas. The terms surviving these two rounds are added by 
a linear combination to separately adjust the right ascension and 
declination offsets by least squares. Table 5 reports the number 
of terms surviving the rounds. Notice the important decreasing 
rate from the first to the second round, which justifies the rounds 
scheme. On the other hand, there is no dropping from the second 
round (one by one regimen) to the third round (least squares fit 
of all validated terms together), what shows that they are truly 
statistically significant. 

Table 6 reports the coefficients derived at the final adjust- 
ment for each family of quasar catalog and stellar frame. As 
expected, the coefficients do not vary from the second round 
to the final adjustment, indicating that no secondary harmonic 
survived. There is only one significant term including the DR5. 
It relates to the 2MASS local correction and depends solely on 
declination (3 rd order). This reflects the high quality of the DR5 
positions and indicates that the local correction did not induce 
further derail upon the input quasar catalogs. Concerning the 



UCACN there is also only one significant term. It relates to the 
GSC2.3 quasar catalog and is of high right ascension order (8 f/l ). 
The lack of UCACN harmonic terms reflects the confinement 
of the sky region where it was applied. The positions originated 
from the B1.0 required 42 harmonic terms (31 for the right as- 
cension offsets), and those from the GSC23 32 harmonic terms 
(20 for the right ascension offsets). Considering the stellar refer- 
ence frames, the UCAC2 required 32 harmonic terms (25 for the 
right ascension offsets), while the 2MASS required 42 harmonic 
terms (28 for the right ascension offsets). In all cases, the right 
ascension offsets required more harmonic terms. In contrast, the 
Ln(sin5) dominated terms in general correspond to larger coeffi- 
cients, which worked to correct the deviations detected at the A4 
equatorial bias analysis. The magnitude terms appear only to the 
l sf order, in 18 B1.0 related terms and 10 GSC23 related terms. 
Thus, magnitude terms are presented as the catalogs authors cau- 
tioned, but they are less important than the zonal terms. Figure 
14 shows the histograms of the harmonic terms coefficients and 
associated errors. The errors resemble well a Poisson distribu- 
tion. The coefficients peak around -6 mas for the right ascension 
offsets, and around +6 mas for the declination offsets. 
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Fig. 11. Equatorial coordinates distribution of the local astrometric corrections derived by the six parameter model using the 2MASS 
stellar reference frame onto the GSC2.3 quasar positions. The clear points correspond to the quasars for which there are ERF radio 
positions. 



5.2. Local inhomogeneities 

Figures 15 and 16 (open symbols) show that after the correction 
by harmonic series, the optical minus radio residuals are equal to 
zero within the statistical significance of 3 x standard error, that 
is without being attached to the scatter of the optical position 
determinations. This is shown for the pairs combining the B1.0 
(Fig. 15) and GSC (Fig. 16) quasar input lists, and the UCAC2 
(squares) and 2MASS (circles) local stellar frames, which are 
the largest and have the most homogenous sky coverage of the 
combinations used here. 

Notwithstanding this, some clumping is seen, without a clear 
pattern over the sky, that can be dealt with further. The ERF sky 
density is such that the mean distance from any given quasar to 
the closest ERF VLBI position is 1.7°. And, on average, there are 
10 ERF VLBI positions within a radius of 4.4° around any given 
quasar, the largest distance being on average 6.3°. Under such 
conditions for each of the combinations of quasar list and stellar 
reference frame, hence for every quasar a correction for the lo- 
cal inhomogeneities towards the VLBI position can be obtained 
from the average of 10 ERF quasars. Within the ERF quasars ra- 



dius, objects are removed one by one if both their right ascension 
and declination offsets are larger than 2cr. On average, 1.3 ERF 
quasars are removed, what shows the robustness of the local in- 
homegeneities corrections. This procedure is similar to that of 
overlapping circles (Taff et al., 1990), but no overlap is actually 
required since the autocorrelation between adjacent corrections 
fades into a characteristic value at a distance of 20°. Figure 17 
illustrates the weakness of the autocorrelations by taking each 
quasar in turn as a pole source, and plotting the standard devi- 
ation in the inhomogeinity correction within growing distance 
rings. For nearby sources, i.e., very small rings, the inhomo- 
geneities are much the same, and the standard deviation is ac- 
cordingly small. Rapidly though, before reaching the distance of 
20°, the standard deviation already reaches the ceiling of scatter 
characteristic of each family of input catalog and stellar refer- 
ence frame solutions. 

By determining in this way the correction for the local inho- 
mogeneities, it is seen in Figs. 15 and 16 (filled symbols) that the 
clumps are clearly minimized. After applying the corrections by 
harmonic functions the optical minus radio offsets are typically 
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Fig. 13. Declination dependency of the equatorial bias (A4). The plotted averages correspond to bins 10° wide in declination. The 
shaded codes indicate the number of ERF quasars within the bins. 
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Fig. 14. Distribution of the coefficients and associated errors of 
the harmonic functions fitting the Aacosd and the A6 residuals. 
All significant terms (first round - see text) are counted here. 



null within 1.3cr, what is reduced further to within 0.7 cr after 
applying the complementary local inhomogeneities corrections . 

6. The LQRF catalog 

The LQRF catalog right ascension and declination equatorial co- 
ordinates are calculated from the values obtained in the final step 
of the data treatment (removal of zonal bias) by a weighted aver- 
age. In the average, for each quasar, the solutions from the three 
input lists (B1.0, GSC23, and DR5) and the three stellar refer- 
ence frames (UCAC2, UCACN, and 2MASS) are combined. The 
weights are the inverse of the square root sum of the internal and 
external errors at each of the correction steps plus the formal er- 
rors of each entry quasar input list. In order though to not carry 
too high or too low a weight for any given combination, for each 
of them separately the errors have been assigned by quartiles. 
It is worth noticing that solutions derived with different weight- 
ing schemes or even without any weighting agree within 10 mas. 
This indicates that the goal of homogenizing the different inputs 
was accomplished. 

Since the final positions originate from the average of multi- 
ple combinations, attaching the final error of them to the squared 
sum of the contributing errors would in fact privilege the sources 
with fewer solutions obtained. Likewise, if attaching the error to 
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Table 4. Direction cosines relative to the ERF after the application of the local astrometric corrections 11 . 



Solution 


N 


Al 


A2 


A3 


A4 


A5 


BU6 


2616 


-33.7 + 6.3 


-8.0 + 6.1 


50.8 ±5.3 


40.5 ± 10.7 


2.6 ±5.3 


BN6 


244 


-5.5 ± 12.9 


31.8 ± 14.1 


123.0 ± 25.2 


1.5 ± 12.8 




BT6 


3764 


-23.8 ±5.1 


1.5+4.7 


23.7 ± 4.6 


27.5 ± 9.7 


-21.0 ±4.3 


BU4 


2522 


-37.5 ± 6.3 


-4.9 ± 6.2 


48.6 ± 5.4 


42.1 ± 10.7 


0.0 ± 5.4 


BN4 


240 


-0.4 ± 12.9 


34.7 ± 14.1 


118.9 ± 25.4 


-0.3 ± 12.8 




BT4 


3762 


-24.6 ± 5.1 


1.4 + 4.7 


24.3 ± 4.6 


24.7 ± 9.6 


-21.7 ±4.3 


GU6 


3238 


-7.0 ± 5.8 


-3.1 ±5.3 


30.4 ± 4.7 


100.8 ± 9.3 


15.2 ±4.8 


GN6 


300 


-7.2 ± 13.2 


16.3 ± 14.2 


-43.4 ± 25.1 


-4.7 ± 13.1 




GT6 


3940 


-4.0 ± 4.7 


-4.9 ± 4.4 


8.7 ±4.3 


30.2 ± 8.7 


-6.1 ± 4.0 


GU4 


3240 


-7.7 ± 5.8 


-2.5 ± 5.3 


29.4 ± 4.7 


99.1 ±9.3 


14.7 ± 4.8 


GN4 


308 


-7.1 ± 13.1 


15.7 ± 14.0 


-40.4 ± 24.7 


-6.8 ± 12.9 




GT4 


3938 


-3.7 ± 4.7 


-6.0 ± 4.4 


8.7 ± 4.3 


30.1 ±8.7 


-6.3 ± 4.0 


SU6 


480 


-5.3+4.7 


-14.6 ±4.8 


-4.8 ± 3.7 


16.6 ±4.1 




SN6 


74 


-5.0 ± 11.6 


14.8 ± 13.9 


32.0 ± 17.8 


9.8 ± 13.1 




ST6 


600 


-13.9 + 5.5 


-14.6 ± 6.0 


-5.1 ±4.9 


20.1 ±5.3 




SU4 


478 


-5.5 ± 4.5 


-12.8 ±4.6 


-5.2 ± 3.5 


16.3 ± 4.0 




SN4 


74 


-6.2+ 11.4 


13.1 ± 13.8 


33.2 ± 17.6 


9.5 ± 13.0 




ST4 


600 


-14.3 + 5.5 


-14.2 ± 6.0 


-4.7 ± 4.9 


19.2 ±5.3 





a Direction cosines Al, A2, and A3 from the comparison of the equatorial positions as found for each family of solutions and as given by 
the radio interferometry for the ERF objects. The equatorial bias A4 (north of declination -20°) and A5 (south of declination -20°), when 
applicable, are also presented. The solutions nomenclature is formed by one letter designating the quasar position input catalog, followed by 
one letter designating the stellar reference frame, and by the codes 6 and 4 informing whether the algorithm is of the six or four parameter 
models. N is the number of equations used to derive the direction cosines and equatorial biases. All angular values are in milli-arcsec. 
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Fig. 17. Fast fading of the autocorrelation of the local inhomo- 
geneities correction. In the figure is plotted the standard devia- 
tion of the corrections in samples taken in rings of increasing ra- 
dius. The upper lines correspond to the USNO B1.0 and GSC2.3 
families of solutions, while the lower lines correspond the SDSS 
DR5 families of solutions. 

the dispersion in the multiple contributing solutions. It is more 
realistic to derive the final error in the equatorial coordinates 
of each source from the external comparison of the LQRF po- 
sition with the radio position. Since in the construction of the 
LQRF the three correction steps focused on the achievement of 
an homogeneous frame, this is done by taking the average of 
the LQRF minus radio positions for the 10 closest ERF sources 



around each quasar. As seen before, this implies a typical ra- 
dius of 4.4° around any LQRF quasar. Figure 18 presents the 
histograms of the errors in the coordinates. The distributions can 
be reckoned to Poisson distributions peaking at 139 mas in right 
ascension and at 130 mas in declination. 

The LQRF contains 100,165 sources. On average, every 
source has a neighbor within 10 arcmin. Figure 19 shows the 
sky distribution counts within squares of 10°. Empty boxes are 
found only in the Galactic plane towards the southern hemi- 
sphere. Figure 20 presents a vector map of the offsets versus the 
VLBI positions, in bins of side 30° and at least 5 quasars. The 
mean offset is 32.7 mas and their distribution appears as ran- 
dom, with the largest offsets corresponding to bins of only 5 or 
6 quasars. 

Figures 21 to 24 show the optical minus radio coordinates 
offsets. The comparison with the equivalent Figs. 5-7 illustrates 
clearly the evenness of the ICRF representation produced by the 
LQRF positions. 

The LQRF catalog is available through 
CDS access (in electronic form via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/VizieR?-source=I/313). 
An extract of its first page appears in Fig. 25. Its columns 
contain 

1 . the J2000 right ascension (h,m,s); 

2. the J2000 declination (d,',"); 

3. the external right ascension error (x cos<5) (mas); 

4. the external declination error (mas); 

5. the number of radio interferometry neighbor positions used 
to evaluate the external errors; 

6. the R magnitude, taken from the LQAC; 

7. the redshift, taken from the LQAC (blank if not defined); 
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Table 5. Number of significant harmonic terms in the three rounds of testing 2 . 
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Harmonic terms surviving in the three rounds of testing for the most significant ones. The three rounds are explained in the text, and the 
third round corresponds to the adopted harmonic fit. All procedure was done independently for the right ascension and declination residuals 
between the solutions and the ERF radio interferometry positions. The solutions nomenclature is formed by one letter designating the quasar 
position input catalog, followed by one letter designating the stellar reference frame, and by the codes 6 and 4 informing whether the algorithm 
is of the six or four parameter model. 



8. the right ascension offsets to the radio interferometry posi- 
tion (blank if there is no radio position) (x costf) (mas); 

9. the declination offsets to the radio interferometry position 
(blank if there is no radio position) (mas); 

10. the J2000 radio interferometry right ascension (h,m,s) (blank 
if there is no radio position); 

11. the J2000 radio interferometry declination (h,m,s) (blank if 
there is no radio position). 

7. Summary and perspectives 

The LQRF is build to derive an optical representation of the 
ICRF, retaining the ICRS concept of being defined by extra- 
galactic objects. Starting from the LQAC entries, quasars were 
indentified in the USNO B1.0, GSC2.3, and SDSS DR5 cata- 
logs. The positions there assigned were made homogeneous with 
respect to the HCRF by applying local corrections obtained in 
small neighborhods around the quasars, relative to the UCAC 
and 2MASS catalogs. The global orientation was then adjusted 
to the ICRF and the zonal departures were removed. The ICRF 
was represented by the ERF containing accurate long base in- 
terferometry radio positions, collected from the ICRF-Ext2, the 
VCS6, and the VLACalib, and selected to precision better than 
10 mas. 

The final LQRF J2000 equatorial cordinates were derived by 
weighted averages of the input catalog positions, once locally, 
zonally and globally corrected as indicated. The errors assigned 
to the LQRF coordinates reflect the departure from the ERF . 

The Large Quasar Reference Frame (LQRF) formed in this 
way contains 100,165 objects, of which 2,142 have accurate 
radio interferometric positions. The comparison between the 
LQRF positions and the ERF positions indicates that the overall 



orientation towards the ICRF, as represented by the Euclidean 
direction cosines is Al=+2.1 + 3.8 mas A2=-0.9 + 3.5 mas 
A3=-2.6 + 3.4 mas with zero equatorial bias to the level of 
2.9 mas . The average offsets to the ICRF are Aacosd=+2.7 
+ 2.9 mas and A<5=+0.3 ± 2.9 mas , with standard deviations 
of <r a = 134.3 mas and cr s - 13 1.1 mas . The internal errors are 
well described by a Poisson representation, peaking at 139 mas 
in right ascension and at 130 mas in declination.. The LQRF 
is planned to be maintained and enlarged in the future, as new 
versions of the LQAC and other quasar input catalogs appear, 
as well as to incorporate newer versions of the UCAC2 stellar 
reference frame. 
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Fig. 25. Example extract of the LQRF catalog. Its columns contain the elements as follows. The J2000 right ascension (h,m,s). The 
J2000 declination (<!,',")■ The external right ascension error (x cos5) (mas). The external declination error (mas). The number of 
radio interferometry neighbor positions used to evaluate the external errors. The R magnitude, taken from the LQAC. The redshift, 
taken from the LQAC (blank if not defined) The right ascension offsets to the radio interferometry position (blank if there is no radio 
position) (x cose)) (mas). The declination offsets to the radio interferometry position (blank if there is no radio position) (mas). The 
J2000 radio interferometry right ascension (h,m,s) (blank if there is no radio position). The J2000 radio interferometry declination 
(h,m,s) (blank if there is no radio position). 
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Fig. 19. Sky density of the LQRF. The counts are in bins of 10° . 
The void regions shown by white circles are in the galactic plane. 
The densest patches lie on the SDSS region. 



Fig. 18. Distribution of the external errors assigned to the right 
ascension and declination J2000 positions of the LQRF sources. 
The distributions can be reckoned to Poisson distributions peak- 
ing at 139 mas on right ascension and at 130 mas on declination. 
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Fig. 20. Vectorial distribution of the systematic deviations (north 
up, east right) of the LQRF to the ICRF within bins of 10° . The 
average value is 32.7 mas, graphically represented by the off map 
arrow. All large values however befall to bins with 5 or 6 sources 
only. 



Fig. 22. Distribution on right ascension of the Ad LQRF minus 
ICRF residuals. There is no systematics above 4 mas and the 
median is 105.3 mas. 
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Fig. 21. Distribution on right ascension of the AacoscJ LQRF mi- 
nus ICRF residuals. There is no systematics above 4 mas and the 
median is 108.3 mas. 
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Fig. 23. Distribution on declination of the Aacostf LQRF minus 
ICRF residuals. There is no systematics above 4 mas and the 
median is 108.3 mas. 
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Fig. 24. Distribution on declination of the A<5 LQRF minus ICRF 
residuals. There is no systematics above 4 mas and the median 
is 105.3 mas. 
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Table 6. Values of the coefficients of the adjusted harmonic functions 3 . . 
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-23.0 ± 6.1 


GU6 


A<5 





2 


A 

4 


+ 1 


-18.6 ±5.7 


GU6 


AS 


1 











+ 7.3 + 2.3 
















GN6 


A<5 








9 


-l 


-33.9+1,1 






n 


n 
u 


1 
1 


+ 1 


4-99 4 + 40 


GT6 


Accost 





1 








+26.0 ± 4.6 


GT6 


AacosS 





3 








-21.0 ± 4.4 


GT6 


A rvc ncA 


I 


I 


1 
1 


— 1 


- 7.9 + 2.2 


GT6 


AS 


1 





u 


u 


+ 6.8 + 2.1 


GT6 


AS 


1 


2 


4 


+1 


+ 7.0 + 2.1 


DU6 


AacosS 












DU6 


AS 












DN6 


Accost 












DN6 


AS 












DT6 


AacosS 





4 








-14.8 + 3.8 


DT6 


AS 













As the results from the six constants and four constants local astrometreic solutions are quite alike, only the results from the six constants 
solutions are displayed. The code for the solutions are as in the previous tables and the coefficients are in milli-arcsec. 



